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The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state is a novel superconducting state in 
a strong magnetic field characterized by the formation of Cooper pairs with nonzero total 
momentum (fc ], —k-\-q |), instead of the ordinary BCS pairs (fc |, — fc J,). A fascinating aspect 
of the FFLO state is that it exhibits inhomogeneous superconducting phases with a spatially 
oscillating order parameter and spin polarization. The FFLO state has been of interest in 
various research fields, not only in superconductors in solid state physics, but also in neutral 
Fermion superfiuid of ultracold atomic gases and in color superconductivity in high energy 
physics. In spite of extensive studies of various superconductors, there has been no undisputed 
experimental verification of the FFLO state, mainly because of the very stringent conditions 
required of the superconducting materials. Among several classes of materials, certain heavy 
fermion and organic superconductors are believed to provide conditions that are favorable to 
the formation of the FFLO state. This review presents recent experimental and theoretical 
developments of the FFLO state mainly in heavy fermion superconductors. In particular we 
address the recently discovered quasi- two-dimensional superconductor CeCoIus, which is a 
strong candidate for the formation of the FFLO state. 

KEYWORDS: FFLO state, upper critical field, vortex state, heavy fermion superconductor, higher Landau 
level, phase diagram 



1. Introduction 

Unconventional superconductivity with a nontrivial 
Cooper pairing state has been the subject of great the- 
oretical and experimental interest. Usually, unconven- 
tional superconductivity is characterized by anisotropic 
superconducting gap functions, which often have zeros 
(nodes) on the Fermi surface. Realization of these nodal 
gap functions has been proved or strongly suggested 
for several classes of materials, including heavy-fermion, 
high-Tc cuprates, and organic superconductors.^"'' These 
unconventional superconductors are characterized by a 
superconducting order parameter which changes sign in 
momentum space along certain directions in the Brillouin 
zone. 

A novel superconductivity that belongs to another 
class of unconventional superconducting state was pre- 
dicted by Fulde and Ferrell,^ and Larkin and Ovchin- 
nikov^ (FFLO) some 4 decades ago. In the FFLO state, 
the superconducting order parameter changes sign in 
real space. The FFLO state originates from the para- 
magnetism of conduction electrons. In the FFLO state, 
the Cooper pair formation occurs between Zeeman split- 
ted parts of the Fermi surface and a new pairing state 
{k t, —k-\-q I) with finite q is realized; the Cooper pairs 
have finite center-of-mass momenta. This contrasts with 
(fc "f , — fc I) pairing in the traditional BCS pairing state. 
Because of the finite q = \q\, the superconducting or- 
der parameter A(r) cx (?/;|(r)-0|(T")) has an oscillating 
component e"' '' and a textured superconducting order 
parameter with a wave length of the order of the super- 
conducting coherence length ^ develops. In the FFLO 
state, spatial symmetry breaking originating from the 
appearance of the q- vector appears, in addition to gauge 



symmetry breaking. As a result, an inhomogeneous su- 
perconducting state appears. Several fascinating super- 
conducting properties which have never been observed 
for the BCS pairing state have been predicted for the 
FFLO state.'^ 

In spite of the straightforward nature of the theoreti- 
cal prediction, this long sought for superconducting state 
has only been reported experimentally recently. This is 
mainly because its occurance requires several very strin- 
gent conditions on the superconducting materials. In 
most type-II superconductors, destruction of supercon- 
ductivity occurs through the orbital pair breaking effect, 
leading to the emergence of the Abrikosov vortex state. 
However, such an orbital effect is always detrimental to 
the FFLO state. ^ Hence, for the FFLO state to occur, the 
orbital pair breaking effect must be weak relative to the 
Pauli paramagnetic effect, so that superconductivity sur- 
vives up to the Pauli limit. Moreover the system needs to 
be very clean, since the FFLO state is readily destroyed 
by impurities.^' Very few superconductors fulfill these 
necessary conditions for the FFLO state. Some of these 
conditions are satisfied in heavy fermion superconduc- 
tors with very large orbital limiting fields and quasi-two- 
dimensional organic superconductors in a magnetic field 
applied nearly parallel to the conducting layers. ^^"^^ In 
addition, anisotropics of the Fermi surface and super- 
conducting gap function are favorable for stability of the 
FFLO state. 

This review presents recent experimental and theoreti- 
cal developments of the FFLO state, addressing the issue 
of its realization in heavy fermion superconductors. In 
particular, we focus on the recently discovered quasi-two- 
dimensional superconductor CeCoIns, for which there is 
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growing experimental evidence of an FFLO state. 

2. FFLO state 

Generally, in type-II superconductors, an applied mag- 
netic field destroys the superconductivity in two dis- 
tinct ways, orbital and paramagnetic pair-breaking ef- 
fects. The orbital effect leads to the emergence of an 
Abrikosov vortex state, forming a regular array of vor- 
tex lines parallel to the field. Then, the kinetic energy 
of the superconducting currents around the vortex cores 
reduces the superconducting condensation energy. The 
orbital limiting field is simply given as H°2^ = $o/27r^^ 
at which vortex cores begin to overlap. Here $0 = Trhc/\e\ 
is the flux quantum. Experimentally, H°2^ at T = is 
commonly derived from the slope of the determined H-T 
phase boundary at as^® 



H°^'°{T = 0) = -0.7 [K] 



dH, 



c2 



dT 



(1) 



The paramagnetic pair-breaking effect originates from 
the Zeeman splitting of single electron energy levels. 
When a magnetic field is applied in the normal state, 
electrons are polarized because the Fermi surface is split 
into spin-up and spin-down electron parts due to the Zee- 
man effect (Pauli paramagnetism). In contrast, in the su- 
perconducting state, spin-singlet Cooper pairs with van- 
ishing spin susceptibility are not spin polarized. Then, 
in order to polarize condensed electrons, the pairs must 
be broken. This destruction of the superconductivity by 
Pauli paramagnetic pair breaking occurs when the Pauli 
paramagnetic energy Ep = ^XnH'^ coincides with the 
superconducting condensation energy Ec = ^A^(0)A^. 
Here Xn = 5(^/^8)^-^(0) is the spin susceptibility in the 
normal state, where g is the spectroscopic splitting factor 
of an electron {g = 2 for a free electron), /xb is the Bohr 
magneton, A is the superconducting gap and A''(0) is the 
density of states at the Fermi level. The Pauli limiting 
upper critical field is estimated to be i?^ = ^A/g/iB 
(Chandrasekhar-Clogston limit) 

Generally, the upper critical field is determined by 
both the orbital and paramagnetic effects. The relative 
importance of the orbital and paramagnetic effects in the 
suppression of the superconductivity is described by the 
Maki parameter^^ 



rrorb 



(2) 



which is the ratio of if"!^ ^c2 zero tempera- 
ture. Since a is of the order of A/ep , where ep is the 
Fermi energy, a is usually much less than unity, indi- 
cating that the influence of the paramagnetic effect is 
negligibly small in most superconductors. However, in 
materials with heavy electron effective mass, in which 
the Fermi energy is small, or in layered materials in a 
magnetic field parallel to the layers, a can be even larger 
than unity. 

2.1 FFLO state in the pure Pauli limit 

We first discuss the superconducting phase of Pauli 
paramagnetically limited superconductors, in which the 
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Fig. 1. H—T phase diagram of 2D and 3D superconductors in 
the pure Pauli limit with no orbital effect. <; = 2 is assumed. 
For a 2D superconductor, H is applied parallel to the plane. 
Below rt (= O.SeTc), the transition becomes first order.^^'^^ 
The broken line represents the metastable transition line. 



uniform magnetic field acts only on the spins of the con- 
duction electrons and all orbital effects are neglected, 
corresponding to the limiting case of infinitely large Maki 
parameter. It has been shown that the Pauli paramag- 
netic effect leads to some peculiar effects on the H-T 
phase diagram. 

For simplicity, we discuss isotropic 2D and 3D super- 
conductors with cylindrical and spherical Fermi surfaces, 
respectively, assuming s-wave superconducting gap sym- 
metry. A 2D superconductor can be realized in a thin 
film in a magnetic field applied parallel to the plane of 
the superconductor. When the film thickness is smaller 
than the coherence length, the orbital pair-breaking ef- 
fect is absent. 

Figure 1 displays the H-T phase diagram of 2D and 
3D superconductors in the Pauli limit with no orbital 
effect. ^^'^^ The initial slope of i/^ at Tc is infinite and 
dHc2/dT cx v/Tc/(Tc - T) in the vicinity of Tc, in con- 
trast to the orbital limited case with finite slope. For 
the orbital limiting case with no Pauli effect, the tran- 
sition at H°^ is of second order. In contrast, for the 
Pauli limiting case, the upper critical field, assuming a 
second order transition line obtained by solving the gap 
equation, decreases with decreasing temperature below 
Tt ~ O.SeTc; i.e. the second order line at T < repre- 
sents a metastable state. Then, below T'^ , the transition 
becomes first order at the upper critical field down to 
T = 0. The upper critical field is H^^. = V^l^/gnB- 

FFLO pointed out that the critical field can be further 
enhanced; the normal state is unstable with respect to 
the second order transition to an inhomogeneous super- 
conducting state. The Pauli paramagnetic pair-breaking 
effect is redticed by the formation of a new pairing state 
{k ],—k + q i) with q ~ 2iJ,BH/hvF, where vp is the 
Fermi velocity, between the Zeeman splitted parts of the 
Fermi surface; spin-up and spin-down electrons can only 
stay bound if the Cooper pairs have finite center-of-mass 
momenta. This contrasts with (fe |, — fc |)-pairing in the 
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Fig. 2. Schematic figure of pairing states, (a) BCS pairing state 
(fc 1,-k + q I), (b) FFLO pairing state (fc T,-k + q I). The 
inner and outer circles represent the Fermi surface of the spin 
down and up bands, respectively. The electron with — fc' + q | is 
not on the inner Fermi surface. 



BCS state. Figure 2 illustrates these pairing states. Here- 
after, we denote the BCS pairing state with q = as the 
BCS state and the pairing state with q 7^ as the FFLO 
state. 

Compared to the BCS state, the FFLO state has a 
smaller condensation energy, butreduction of the Zee- 
man energy can stabilize the FFLO state. Then upper 
critical field in the FFLO state i/pFLO can be higher 
than H^2- At T = 0, H^^^q ~ f.5lAo/gAiB ^ 1-07 H^^ 
for 3D,5'6,29,32 iJl^LO - SAo/g/iB ^ l-42i7P2 for 
13,14,33,34 Pqj. -^j^ TJi^LO diverges as T -> 0, as wiU 
be discussed later. '^^^^^ It turns out that the FFLO state 
only appears at T < T^; the point (T^, H^) is a tricritical 
point, at which the normal, uniform BCS superconduct- 
ing, and FFLO phases all meet. 

The q-vector in the FFLO state gives rise to spatial 
symmetry breaking. As a result, the superconducting or- 
der parameter exhibits spatial oscillation. A possible su- 
perconducting order parameter, originally proposed by 
Fulde and Ferrell, is 

A(r) = Ale"^■^ (3) 

in which the amplitude of the superconducting order pa- 
rameter is homogeneous, but the phase changes in real 
space. In this state, depairing occurs for some part of 
the Fermi surface (the shaded portion in Fig. 3) and 
pairing can ocur in the remaining regions of the Fermi 
surface.^' ■^^ 

In the gap equation, the two solutions Aic"' '^ and 
Aie^"^ '' represent degenerate superconducting states at 
a transition for a given q. The degeneracy is lifted by 
forming the linear combinations of e'^ '' and e""' '". It 
has been found that 

A(r) = Ai(e'«-'^ + e""'-'^) = 2Ai cos(q • r) (4) 

provides a lower free energy than eq. (3) for the whole 
temperature range below T^.^ This state was originally 
proposed by Larkin and Ovchinnikov. In this state, the 
superconductor becomes spatially inhomogeneous, be- 
cause the order parameter undergoes ID spatial oscil- 
lation with a wave length of 27r/g and, as a result, the 
normal quasiparticle and paramagnetic moments appear 
periodically. 

Generally, because of the symmetry of the system, 
there are more than two equivalent q-vectors which give 
the same upper critical field. Then, the order parameter 




Fig. 3. Pairing state described by eq. (3). The shaded region is 
the blocking region where Cooper pairs are not formed. 



is expressed by a linear combination of e"''" '^, 

M 

A(r) = ^rne"'"''^, (5) 

m— 1 

where M denotes the number of equivalent q- 
vectors.®''^^"''^ In an isotropic system, M = oo. 

Figure 4 displays the H-T phase diagram and quasi- 
particle structure in a 2D thin film in a parallel field. 
Just below T\ the order parameter described by eq. (4) 
is stabilized. However, as the temperature decreases, or- 
der parameters with linear combinations of larger num- 
ber of plane waves stabilize. '^^ This is because the in- 
creasing number of plane waves increases the number of 
nodes, which in turn reduces the Pauli paramagnetic en- 
ergy of the quasiparticles excited around the nodes. At 
lower temperatures, the "triangular state" , with a gap 
function of the form 

A(r) = Ai (c'"!-'- + + e'"^ ''), (6) 

has a lower free energy, with q^ = ((7,0,0), q2 = 
i-q/2, \/3g/2, 0), and qg = i-q/2, -V3q/2, 0), where we 
take the z-axis to be perpendicular to the 2D plane. A 
further decrease of the temperature stabilizes the "square 
state" , with a gap function of the form 

A(r) = Ai {cos{qx) + cos{qy)), (7) 

and the "hexagonal state" , with a gap function of the 
form 

A(r) = Ai (cos(qj • r) + cos(q2 • r) + cos{q^ ■ r)). (8) 

In a 3D isotropic system, a recent calculation**^ demon- 
strates that the first order transition line is slightly higher 
than the second order transition line at all temperatures 
below . Moreover, similar to the 2D case, the order pa- 
rameter is cos(q • r) just below T\ but at lower tempera- 
tures the order parameter switches to more complicated 
structures which are expressed by a sum of two or three 
cosines. It has also been suggested that at zero tempera- 
ture a state, which is expressed by a sum of eight plane 
waves forming a face-centered cubic structure, appears.^* 

We note that the FFLO state had not been reported 
for thin films until now. The main reason for this is that 
the FFLO state is readily destroyed by impurities,^' 
as will be discussed later. In fact, most thin films are in 
the dirty limit, £, ^ £, where i is the mean free path of 
the electrons. Even in clean films, surface scattering may 
destroy the FFLO state. Even if a very clean thin film 
could be produced, the magnetic field must be applied 
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Fig. 4. H—T phase diagram of a 2D superconductor in a parallel 
field. The color regions are in the FFLO state. (Tt.fft) is the 
tricritical point. The solid and dashed lines represent the second 
and first order transition lines, respectively. Just below T^, the 
order parameter undergoes ID spatial oscillation, as expressed by 
eq. (4). At lower temperatures, triangular, square, and hexagonal 
FFLO states, expressed by eqs. (6), (7) and (8), respectively, 
appear. For details, see the text. The quasiparticle densities of 
the states in each FFLO state are shown. 



very precisely parallel to the plane to produce the FFLO 
state. Otherwise a small but finite perpendicular field 
component leads to the formation of a vortex state, which 
destroys the FFLO state (see §2.2). 

The stability of the FFLO state with respect to ther- 
mal fluctuations is different from that of BCS supercon- 
ductors.^^' ^■^ In an isotropic system, the direction of the 
q-vector is arbitrary. For 3D systems, long-range-order 
(LRO) is destroyed by thermal fluctuations of the q- 
vector and quasi long range order (QLRO), characterized 
by a correlation function with a power law decay, is es- 
tabhshed.^^'^^ For 2D systems, it is well known that LRO 
is absent even in the BCS state, and superconductivity 
with QLRO is achieved through a Kosterlitz-Thouless 
transition. In the FFLO state with ID oscillation ex- 
pressed by eqs. (3) and (4), even QLRO is absent due to 
directional fluctuation of the q- vector. ' ^'^ However, it 
has been shown that LRO is established even in an FFLO 
state in the presence of crystal anisotropy,^^ which exists 
in real systems, because the anisotropy fixes the q-vector 
in a certain direction. It has also been shown that the 
fluctuations lead to a divergent spin susceptibility and a 
breakdown of the Fermi-liquid behavior at the quantum 
critical point. Moreover, in the FFLO state of quasi-2D 
c?-wave superconductors, a fluctuation-driven first-order 
transition has been proposed. 

2.2 Orbital pair-breaking 

In real bulk type-II superconductors, the orbital effect 
is invariably present. Gruenberg and Gunther examined 
the stability of the FFLO state against the orbital ef- 
fect in 3D isotropic systems with s-wave pairing.* It has 
been shown that orbital effects are detrimental to the 
formation of the FFLO state. The FFLO state can exist 



at finite temperatures if a is larger than 1.8. The FFLO 
region shrinks considerably from that in the absence of 
the orbital effect. 

In the presence of a magnetic field H = (0, 0, iJ), the 
spatial variation of the order parameter is determined 
by the cyclotron motion of the Cooper pairs in the plane 
perpendicular to H . Generally, the orbital effect forces 
the solutions for the order parameter to be eigenfunctions 
of the operator iT^, where U — — iftV — with the 
vector potential A — (0, Hx, 0). This problem is equiva- 
lent to that of a charged particle in a constant magnetic 
field. The operator iT^ can be rewritten as 

2h^ 



^2 



1 Fl'^ 
2> 



(9) 



where rj^ and rj are the boson creation and annihilation 
operators defined by 



V 



t - 



V2n 

Ail 
V2h 



(n. - in,) 
(n, + inj,). 



(10) 



respectively. Here, = \J $o/27r_ff. Eigenfunctions with 
respect to the degrees of freedom perpendicular to H are 
described by the Abrikosov functions with Landau level 
index n, 



0W(p) = (-l)V^-^iJ„[y2 



(11) 



for any real number k and Xk = k^n, where _ff„ is the 
Hermite polynomial and p = {x,y). The lowest Landau 
level solution (n = 0) is given by 

(x-Xk)'^. 



</)},"-^(p)=e"=^exp[- 



(12) 



In the BCS state, the eigenvalues of the operator iT 



2^2 1 



(13) 



Usually, the BCS state is described by the lowest Landau 
level solution 



A(r)(x0f (p), 



(14) 



(15) 



which gives the highest upper critical field. 

The coexistence of both the FFLO state and vortex 
modulation has been investigated through exact calcula- 
tions of the upper critical field assuming that the tran- 
sition is of second order.®' ^^'^"^ In the FFLO state, the 
eigenvalues of the operator iJ^ are 

2^ 

The second term represents the kinetic energy of the 
Cooper pairs along H. We note that the increase of the 
kinetic energy in eq. (15) is compensated by the reduc- 
tion of the Pauli paramagnetic energy. Assuming that the 
lowest Landau level gives the highest upper critical field 
even in the FFLO state, the order parameter is given by 

A(r)cxexp(i(7z)4''(p) (16) 
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Fig. 5. H—T phase diagram in the presence of the orbital effect, 
as proposed in Ref. [8] The transitions from the normal state to 
the FFLO and BCS states are of second order (solid line), while 
the transition from the FFLO state to the BCS state is of first 
order (thick solid line). 



or 

A(r)acos(qz)4"^(p). (17) 

As indicated by eq. (17), the order parameter exhibits 
a modulation with wave vector q parallel to H and 2D 
nodal sheets appear perpendicular to H . Equation (17) 
clearly indicates that parallel motions of the Cooper pairs 
create the 2D nodal sheets. 

The H-T phase diagram in the presence of the orbital 
effect proposed by Gruenberg and Gunther is schemat- 
ically shown in Fig. 5. When the orbital effect is not 
weak, the transitions from the normal state to the FFLO 
and BCS states are of second order, while the transition 
from the FFLO to the BCS state is of first order. Recent 
calculations indicate that the phase diagram is modified 
in some cases. ^^"^^ For example, according to Ref. [75], 
in which the 4th and 6th order terms of the Ginzburg- 
Landau free energy are incorporated, some portion of the 
transition line at the upper critical field is of first order, 
while the transition from the FFLO state to the BCS 
state that branches from the first order line is of second 
order, as shown in Fig. 6. 

A further increase in the relative importance of the 
Pauli paramagnetic effect produces an FFLO state 
with the Abrikosov functions of a higher Landau 
level. This is because the Pauli paramagnetic 
energy is reduced as the number of nodes is increased, 
which in turn increases the number of polarized quasipar- 
ticles in the vicinity of the nodes. This situation would 
be realized in heavy fermion superconductors with very 
large Maki parameter and in 2D layered superconduc- 
tors in magnetic fields slightly tilted from the 2D plane. 
In a 2D system, very small perpendicular component of 
H leads to a very large effective Maki parameter.^'^^^^'^^ 
Figure 7 shows a phase diagram of a 2D system in a 
slightly tilted magnetic field, obtained in Ref. [77] . A sim- 
ilar phase diagram has also been obtained in Ref. [78]. 
A rich variety of quasiparticle structures associated with 
the higher Landau level indexes appear as a function of 



1 

T/T, 

Fig. 6. H—T phase diagram for a 3D system with large Maki pa- 
rameter (a = n denotes the Landau level index. The thick 
and thin solid lines represent the first and second order transition 
lines, respectively. 



the tilted angle and temperature.™"^^ For a 3D system, 
it has been shown that an FFLO state can coexist with 
the Abrikosov state with higher Landau levels. Fig- 
ure 6 schematically shows the H-T phase diagram for 
a 3D system with large Maki parameter as proposed in 
Ref. [75]. At very low temperatures, the Abrikosov func- 
tion with the first Landau level {n — \) provides the 
highest upper critical field. In this state, nodal planes 
are formed periodically parallel to the vortex lattice. 

2.3 Anisotropics oj Fermi surface and gap Junction 

As discussed in the previous sections, the FFLO state 
is restricted in an extremely narrow region in the H- 
T plane (-ffpFLO ~ 1-07 x H^2 tlie absence of the 
orbital pair-breaking effect). In this circumstance, the 
FFLO state is readily destroyed by several effects, such 
as internal field enhancement due to electron correla- 
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T/T, 

Fig. 7. Upper critical fields of states with n of a 2D system in a 
magnetic field slightly tilted from the plane of the superconduc- 
tor.'^'^ n denotes the Landau level index. The transitions at Hc2 
are assumed to be of second order. The physical critical field is 
the maximum one at each temperature. The broken and dotted 
lines show the critical fields of the BCS state (q = 0) and the 
FFLO state (q 0) in a parallel field, respectively. 
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tions (Fermi liquid correction). -"^'^'^ In fact, a very 
slight enhancement of the uniform susceptibility x from 
that without correlations xo, x/xo = 1-04, is large 
enough to suppress the FFLO state. -'^'^-'^ For stability of 
the FFLO state in real systems, it is necessary that the 
Fermi surface or the gap function be anisotropic. These 
anisotropies stabilize the FFLO state through optimiza- 
tion of the q- vector. 

The stability of the FFLO state also crucially depends 
on the dimensionality of the system. This can be under- 
stood from Fig. 2. The FFLO state is more stable when 
a large area of the spin-up Fermi surface is connected to 
the spin-down surface by the q-vector. This bears some 
analogy with the nesting effect of charge- density- wave 
(CDW) and spin-density-wave (SDW) transitions.^^' 
For a 3D system with a spherical Fermi surface, the spin- 
down and spin-down Fermi surfaces touch only at a point 
by a shift of the g- vector. For a 2D system with a cylindri- 
cal Fermi surface, the two Fermi surfaces touch on a line 
by a shift of the g-vector. A particular case is a ID sys- 
tem with a flat Fermi surface where the whole region of 
the Fermi surfaces touch by a single q- vector. Therefore, 
low dimensionality favors the FFLO state, extending the 
FFLO regime. However, for the ID case, CDW or SDW 
transitions would occur at a finite temperature and the 
Fermi surface disappears. Therefore, the best system to 
observe the FFLO state appears to be a 2D system. 
Even in a 3D system, if the Fermi surface has sonic; flat 
regions, the FFLO state is expected to be stabilized when 
H is oriented in an optimal direction of the g-vector. 

The above argument considers the FFLO state for 
superconductors with isotropic s-wave symmetry. Most 
of the candidate compounds with large Maki parame- 
ter belong to two classes of materials, quasi-2D organic 
and heavy fermion superconductors. All these materi- 
als have been argued to be unconventional superconduc- 
tors with anisotropic gap fimctions. most likely to have 
nodes, and in this way differ from the case originally con- 
sidered by FFLO. Hence, the experimental identification 
of nodal structures has stimulated theoretical interest in 
the FFLO state in anisotropic d-wave pairing symme- 

^j.y 43-45, 74, 77, 102 

It has been shown that the stability of the FFLO state 
also depends on gap anisotropy. The gap function de- 
scribing the FFLO state of anisotropic pairing has the 
form, for example, 

oo 

A(r, fe) (X 7(fc) cos{qz) ^ C„0(f ' (p), (18) 

with k = k/k, where r = (p, z) is the ccnter-of-mass co- 
ordinate and k is the relative momentum of the Cooper 
pair.''*'''^ Here, 7(fe) describes the momentum depen- 
dence of anisotropic pairing. For d,j.2_y2-wme pairing, 
7(fc) oc A:^ — ky, and for dxy, 7(fe) oc k^ky. For d-wave 
symmetry in a 2D system with a cylindrical Fermi sur- 
face, the FFLO state is more stable when the g-vector 
is oriented in the direction of maximum gap at low 
temperatures. In actual systems, the optimal q-vector 
depends on both the Fermi surface and gap function 
anisotropies.''6'i°2 in k-(BEDT-TTF)2X,i°2 the FFLO 



state has been studied by taking into account both ef- 
fects. 

2.4 Impurity Effect 

Impurity effects on the FFLO state have been stud- 
ied by several authors.^' ^'^'^^'^^ It has been shown that 
the FFLO state is unstable against nonmagnetic impu- 
rities, in contrast to the BCS state of s-wave pairing. 
This is analogous with anisotropic superconductivity, in 
which the mixing of order parameters with opposite signs 
in momentum space caused by impurity scatterings sup- 
presses the pair amplitude. Similarly, in the FFLO state, 
mixing by impurity scattering in real space suppresses 
the pair amplitude. 

It has also been shown that non-magnetic impurities 
dramatically change the FFLO phase diagram. In the 
presence of impurities, the behavior of the FFLO phases 
in unconventional superconductors is qualitatively differ- 
ent from that in s-wave superconductors.^^ For example, 
it has been suggested that in d-wave superconductors, 
states with order parameters expressed by eqs. (3) and 
(4) both appear in the same H-T plane. 

A possible change of the order of the phase transition 
from normal to the BCS and FFLO states due to im- 
purities has also been suggested. A rich variety of phase 
diagrams has been proposed when both the orbital pair- 
breaking and impurity effects are taken into account. 

2.5 Quasiparticle structure of a vortex in the FFLO 
state 

The quasiparticle structure in the superconducting 
state in a magnetic field has been studied extensively. In 
the BCS state, the quasiparticles are spatially localized 
within vortex cores forming Caroli-de Gennes-Matricon 
bound states with discrete energy levels. The quasi- 
particle structure in the Pauli limit without the orbital 
effect has been discussed in §2.1. (For the 2D case, sec 
Fig. 4.) The quasiparticle structure in the FFLO state in 
the presence of the orbital effect has been addressed by 
several authors. Recently, the low energy quasi- 
particle excitation and spin polarization spectrum in the 
FFLO state for the isotropic 3D case has been exam- 
ined in the framework of the microscopic Bogoliubov-de 
Gennes theory, taking into account the orbital and Pauli 
paramagnetic effects in an equal footing. Figures 8(c) 
and 8(d) display schematically the quasiparticle excita- 
tions in the BCS and FFLO states, respectively. It is re- 
markable in the FFLO state that the quasiparticles are 
not excited in the region where the (vertical) vortex lines 
intersect the (horizontal) planar nodes. 

This peculiar spatial quasiparticle structure is charac- 
teristic of the topology of the FFLO vortex, as schemat- 
ically shown in Fig. 8(e). Generally, quasiparticles are 
spatially localized at the region where the phase of the 
superconducting order parameter is shifted by zLtt;^^ the 
bound states appear as a result of interference of the su- 
perconducting wave function caused by Andreev reflec- 
tion by pair potentials with opposite signs. In Fig. 8(e), 
(i) and (ii) represent the path cutting through the vor- 
tex line with 27r-phase winding (A oc e'^) and the path 
crossing the FFLO nodal plane (A a sin(g • r), where 
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Fig. 8. (a) Schematic view of the flux hne lattice in the BCS 
state, (b) Flux hne lattice in the FFLO state. In the FFLO state, 
planar nodes appear periodically perpendicular to the flux line, 
which leads to a segmentation of the vortices into pieces of a 
length A = 27r/g. (c) Schematic flgure of the quasiparticle exci- 
tation spectrum in the BCS state. The quasiparticles are excited 
within the vortex cores (blue region), (d) Quasiparticle excita- 
tion spectrum in the FFLO state. The quasiparticles here are 
also excited around the FFLO planar nodes. The quasiparticle 
regions within the vortex lines do not spatially overlap with those 
in the FFLO planar nodes, (e) Schematic view of (vertical) the 
vortex line and (horizontal) the FFLO nodal plane. Quasiparti- 
cles experience sign change of A(r) when tracing either path (i) 
or (ii) through the tt shift of the phase in the order parameter, 
while no sign change occurs when tracing a path (iii). For details, 
see the text. 



we take the FFLO nodal plane to be 2; = 0, respectively. 
Quasiparticles tracing either path (i) or (ii) experience 
the sign change of A(r) through the tt shift of the phase 
in the order parameter. The phase of the order param- 
eter around the vortex line is additionally twisted by tt 
when crossing the FFLO sheet. When tracing the path 
(iii) (the dotted line in Fig. 8(e)), the quasiparticles do 
not experience a sign change of the pair potential, be- 
cause of a 7r± TT phase shift. Therefore, bound states are 
not formed in the region where the vortex line intersects 
with the FFLO sheet. 

The peculiar quasiparticle structure in the FFLO state 
is expected to change the elastic properties of the flux 
hne lattice (FLL)}'^''^° Notably, the occurrence of FFLO 
nodal sheets perpendicular to the vortex lattice is ex- 
pected to lead to a softening of the FLL tilt modulus 
C44. This is due to the fact that a deformation of the 
nodal planes entails a deformation of the FLL (sec §5.3). 
The decreased order parameter in the FFLO state re- 
duces the pinning strength that characterizes the FLL 
interaction with crystal defects. In other words, defects 
in the vicinity of the nodal planes no longer interact with 
the FLL, even though pinning of the nodal planes is pos- 
sibile. Defects between the nodal planes exert a smaller 



pinning force (proportional to |Ap). 

3. Candidate compounds for the FFLO state 

Summarizing the requirements for the formation of the 
FFLO state: 

(1) Strongly type-II superconductors with very large 
Ginzburg-Landau parameter k = A/^ ^ 1 and large 
Maki parameter a, such that the upper critical field 
can easily approach the Pauli paramagnetic limit. 

(2) Very clean, £, £, since the FFLO state is readily 
destroyed by impurities. 

(3) Anisotropics of the Fermi-surface and the gap func- 
tion can stabilize the FFLO state. 

The main reason that the FFLO state has never been 
observed in conventional superconductors resides in the 
fact that the requirements for the formation of the FFLO 
state are extremely hard to fulfil. The question of observ- 
ing the FFLO state in unconventional superconductors 
has been addressed recently. Heavy fermion and organic 
superconductors are attractive classes of superconduc- 
tors for forming the FFLO state because of their poten- 
tially large Maki parameters, very large mean free paths, 
and unconventional pairing symmetries.* 

Heavy electron masses lead to low Fermi velocities of 
the quasiparticles and, in turn, enhance H°2^ and the 
Maki parameter. For this reason, heavy fermion super- 
conductors have had considerable attention. In fact, the 
features of the magnetization and the H-T phase dia- 
grams of CcRu2,"''^^ UPd2Al3,"'i6.57 and UBeis were 
taken as possible signatures of the FFLO state. Very re- 
cently, a very strong candidate for exhibiting the FFLO 
state, CeCoIns, has been discovered for which a new su- 
perconducting phase has been confirmed by the thermo- 
dynamic second order transition. 

In layered organic superconductors, when a mag- 
netic field is applied parallel to the layers, the 
orbital effect is strongly suppressed. ^"'^ Indica- 
tions for a possible FFLO state have been reported 
for A-(BETS)2GaCl4,i°'' A-(BETS)2FeCl4,i°^-"° 
and K-(BEDT-TTF)2Cu(NCS)2.^°^'"^""^ For A- 
(BETS)2GaCl4, a kink in the thermal conductivity 
suggests a transition to the FFLO state. For k- 
(BEDT-TTF)2Cu(NCS)2 a similar feature in the 
magnetization has been identified. ^ ■ In spite of 
several observations which appear to support the FFLO 
state, evidence of the thermodynamic phase transition 
between the FFLO and the BCS-pairing state has not 
been reported. A detailed description of a possible 
FFLO state in organic superconductors is beyond the 
scope of this article. For a recent review, see Ref. [114]. 

An FFLO state in layered high-Tc cuprates with strong 
superconducting fluctuations is highly unlikely, because 
the region of the FFLO state is replaced by a vortex 
liquid state below the mean field upper critical field line, 
which is a nonsuperconducting phase. 

Besides heavy fermion and organic compounds, some 
Chevrel phase materias, like PbMogSs might be consid- 
ered candidates for exhibiting the FFLO state, be- 
cause of their very large upper critical field. However, 
the very small mean free path of these materials {£ ^ ^) 
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prevent the formation of the FFLO state. Borocarbide 
superconductor TmNi2B2C have also been suggested to 
be a promising candidate for the FFLO state because of 
a large paramagnetic effect. 

4. Possible FFLO state in heavy fermion super- 
conductors 

We here address the issue of the superconducting 
states in CeRu2, UPd2Al3, UBeia, and CeCoIns. In these 
materials, a possible existence of the FFLO state has 
been suggested by a number of groups. For the 
following discussion, the physical parameters related to 
the FFLO transition in the normal and superconducting 
states of these materials are listed in Table I. 

4-1 CeRu2 and UPd2Aly,: is the peak effect a signature 
of the FFLO state? 

CeRu2 is a superconducting cubic Laves phase com- 
pound with Tc = 6.1 K. Since its discovery in 1958, 
CeRu2 had been believed to be a conventional s-wave 
type-II superconductor. Recent experiments have re- 
vealed that the superconducting gap function is modu- 
lated at the Fermi surface, Amin/Amax ~ 0.4, indicating 
an anisotropic s-wave.^®'^° 

UPd2Al3 is a heavy fermion superconductor dis- 
playing quite unique properties, including the coexis- 
tence of superconductivity and antiferromagnetic order- 
ing with an atomic size local moment (/i = 0.85/xb), 
and strong coupling between the magnetic excitation 
and superconductivity.^^' Several experiments have re- 
vealed that the superconducting gap symmetry is most 

likely to be d-wave with line nodes parallel to the layers, 
A(fc) = Acosfc^c.4'21.64 

These materials have attracted interest because a pro- 
nounced hysteresis behavior has been observed in a nar- 
row field range below Hc2, Hi < H < Hc2, in several dif- 
ferent physical properties, such as the magnetization, ac 
susceptibility, ultrasound velocities, and thermal expan- 
sion. On the other hand, these quantities are reversible in 
a wide field range below Hi down to the lower critical field 
Hci- Particularly, at Hi < H < Hc2, the isothermal de- 
magnetization process in CeRu2 and UPd2Al3 exhibits a 
hysteresis loop with a sharp peak at H^^ak with decreas- 
ing H and a sharp negative peak at the same field with 
increasing H. This peak eff'ect arises from the critical 
state, indicating the occurrence of a peak in the critical 
current density of the system. 

The origin of the peak effect has been discussed in rela- 
tion to a possible FFLO state, since the vortex lines can 
become segmented into coupled pieces of lengths compa- 
rable to ^, and become flexible in a qualitatively similar 
way as a pancake vortex in high-Tp cuprates^^ near the 
upper critical field. These segments could conform more 
easily to pinning centers, thereby enhancing pinning. In 
addition, the transition line (Tpeak, -ffpeak) in CeRu2 and 
UPd2Al3 resembles the FFLO transition line in the H-T 
plane (see Fig. 5, for example). These features are also 
taken as possible signatures of the FFLO state. More- 
over, the large electron mass and anisotropic gap sym- 
metry of UPd2Al3 are favorable for FFLO formation. 



However, subsequent studies have cast doubt on the 
interpretation of the peak effect in these two compounds 
in terms of the FFLO transition. For instance, the Maki 
parameter of CeRu2 is less than unity, presumably be- 
cause CeRu2 is not a heavy fermion material and hence 
the orbital critical field H°^ is not large (see Table II). 
Moreover, the transition line (Tpcak, -f^peak) splits from 
the ifc2-fine at Tpeak/^c w 0.85 for CeRu2 and w 0.80 for 
UPd2Al3. These values are much larger than the theoret- 
ical predictions in Refs. [5,6,8]. A possible explanation for 
the large value of Tpeak/^c in UPd2Al3 is that the actual 
Tc at zero field is depressed by the pair-breaking effect 
due to the exchange field in the antiferromagnetically or- 
dered state. However, as pointed out in Ref. [66], the 
calculation using the parameters obtained for UPd2Al3 
provide a which is much smaller than 0.56Tc. Thus 
the features of the peak effect in CeRu2 and UPd2Al3 
are clearly at odds with the FFLO transition. 

What then is the origin of the peak effect in CeRu2 and 
UPd2Al3? In the last decade, understanding of the peak 
effect has progressed. The peak effect is observed in clean 
superconductors with weak pinning. Among them, the 
very anisotropic high-Tc cuprate Bi2Sr2CaCu208-i-5 and 
2if-NbSe2 have been widely studied due to their weak 
pinning properties. It is widely accepted that in these 
systems the peak effect is ascribed to a disorder-driven 
first-order transition from the ordered Bragg-glass to the 
disordered vortex glass, caused by competition between 
the elastic and pinning energies. ^''"^^ 

An important question is: Is the peak effect in CeRu2 
and UPd2Al3 of the same origin as those in other clean 
superconductors, i.e., a disorder-driven transition generic 
to weakly pinned vortices? According to subsequent 
studies of ac-susceptibility in 2i?-NbSe2 and CeRu2, 
both materials display remarkable similarities in the peak 
effect region, which can be explained by disorder-induced 
fracturing and entanglement^'^ or transitions from collec- 
tive to single- vortex pinning of the vortex lattice. ''^ Based 
on these results, it has been concluded that the FFLO 
state in CeRu2 is not required to explain the peak effect. 

Thus, subsequent research has cast doubt on the in- 
terpretation of the peak effect in CeRu2 and UPd2Al3 in 
terms of the FFLO state. At the present stage of study, it 
is safe to conclude that the FFLO state is very unlikely in 
CeRu2 and UPd2Al3, though the detailed quasiparticle 
structure in the peak effect region should be scrutinized. 

The observation of a possible FFLO state has also been 
argued for UBeis,^^ which exhibits a strong upturn in 
the upper critical field at low temperatures.^^ However, 
recent NMR Knight shift measurements have revealed 
that the parity of UBeis appears to be spin-triplet with 
no Pauli paramagnetic effect on the superconductivity.^^ 

5. Evidence for an FFLO state in CeCoIns 

5.1 CeCoIn^ 

The recently discovered superconductor CeCoIns is a 
very strong candidate for exhibiting the FFLO state. 
CeCoIus is a heavy fermion superconductor with a 
of 2.3 K, the highest among Ce- and U- based heavy 
fermion superconductors.^^ It exhibits a layered struc- 
ture of alternating conducting Celns planes (Celns is 
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Table I. Physical parameters in the normal and superconducting states of CoR,U2, UPd2Al3, and CeCoIns. The orbital limited upper 
critical fields in CeCoIns are determined by eq. (1). 
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a superconductor under pressure) and less conducting 
CoIn2 planes, suggesting a quasi-2D electronic nature. 
This is supported by de Haas-van Alphen studies which 
reveal quasi-cylindrical sheets of the Fermi surface^^ 
and by NMR relaxation rate measurements which re- 
veal the presence of quasi-2D spin fluctuation.^^ Related 
to the quasi-two-dimensionality, the upper critical field 
for II ah, h1^2i is nearly twice as large as that for 
H _L ah, H^. According to several experiments, in- 
cluding studies of the specific heat,*^ thermal conduc- 
tivity,®^"*^ NMR,*5'89'9° penetration depth,^! and An- 
dreev reflection, the superconducting gap symmetry of 
CeCoIn5 is most likely to be d2;2_j^2-wave with line nodes 
perpendicular to the plane.*'^'^'^'^'' 

The H-T phase diagram of CeCoIns exhibits several 
fascinating properties, which have not been observed in 
other superconductors. In a strong magnetic field, the 
superconducting transition is of first order at low tem- 
peratures for both configurations, H \\ ah and H _L ah, 
as shown in a number of experiments, including jumps 
of the thermal conductivity,^^' ®^ magnetization,^^ pen- 
etration depth, ultrasound velocity,^" thermal expan- 
sion,^'^ magnetostriction,^^ and NMR Knight shift, ^^'^'^ 
and a sharp specific heat anomaly. For H J- ah and 
H II ah, the first order transition terminates at ^ 0.3 Tc 
and ~ 0.4 Tc, respectively, above which the transition 
becomes second order up to Tc. 

The first order transition at the upper critical field 
indicates that the Pauli paramagnetic effect dominates 
over the orbital effect in both field directions. Accord- 
ing to Ref. [75] , for observation of the first order transi- 
tion, the system must be very clean. Quantitatively, for 
a — I, the transition becomes first order only when i is 
more than 10 times longer than ^. These requirements 
are well satisfied in CeCoIns. Indeed, the orbital limit- 
ing fields obtained from eq. (1) are 38.6 T and 17.7 T 
for H \\ ah and H _L ah, respectively. Then, assuming 
that the upper critical fields are close to the Pauli lim- 
iting fileds, H^2 - 11-8 T and 4.95 T for H \\ ah and 
H _L ah, respectively, the Maki parameters determined 
by eq.2 are estimated to be a" ~ 4.6 for H \\ ah and 
— 5.0 for H _L ab.^^ Moreover, microwave and ther- 
mal Hall conductivity measurements demonstrate that 
i is more than 100 times longer than ^ in the supercon- 
ducting state, indicating that CeCoIns is in an extremely 
clean regime.*®'^® 

5.2 FFLO phase in CeCoIn^ 

For H II ah, a second anomaly in the heat capacity 
has been observed within the superconducting state; a 
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Fig. 9. Experimentally determined H—T phase diagrams of 
CeCoIns at low temperatures and high field for (a) H A- ab 
and (b) H \\ ab. The colored portions display the FFLO (pink) 
and BCS (blue) regions. The black and green lines represent the 
upper critical fields, which are first order and second order, re- 
spectively. The red lines represent the phase boundary separating 
the FFLO and BCS states. 



second order transition takes place in the vicinity of 77^2 
at low temperatures.^*' This second-order transition 
line (T|*,iJ|*) branches from the first order ijjj line and 
decreases rapidly with decreasing T , indicating the pres- 
ence of a new superconducting phase. The second or- 
der anomaly {T*^,H*^) has been confirmed by the sub- 
sequent studies of penetration depth, thermal conduc- 
tivity,^^ ultrasound velocity,^" magnetization,^^ magne- 
tostriction^^ and NMR.^**'^* This novel high field super- 
conducting state at the low-T/high-iJ corner of the H-T 
plane has been identified as the FFLO state. 

For H _L ah, heat capacity measurements have also 
identified an anomaly just below H^}^ However, the 
study of the superconducting state has focused on the 
case of II ah. Recently, a new superconducting phase 
for H J- ah has been confirmed by NMR^*^ and magneti- 
zation^^ measurements. From the analogy of the results 
for H II ab, the presence of FFLO state for H J- ab has 
been suggested. 

CeCoIns appears to meet all the strict requirements 
placed on the existence of the FFLO state discussed in 
the beginning of §3. Figures 9(a) and (b) display the H- 
T phase diagrams for H \\ ah and for H _L ah, respec- 
tively, determined from several experiments. The tran- 
sition lines determined by different techniques are very 
close to each other. We next discuss the FFLO state as 
probed by various methods. 
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5.3 FFLO state in a parallel field 

5.3.1 Penetration depth 

The penetration depth in the a6-plane for H \\ ah 
has been measured at low temperatures near Hc2.^^ It 
has been shown that the penetration depth increases at 
(Tipifip. This can be interpreted as representing a de- 
crease of the superfluid density due to the formation of 
the FFLO state. 

5.3.2 Magnetization and magnetostriction 
Magnetization M{H) as a function of H has been 

measured by two groups near _ff|'2.^^'^^ According to 
Ref. [27], a change in the monotonic variation of the 
magnetization M{H) is observed at H'^y This anomaly 
is clearly seen by plotting the field derivative of the mag- 
netization AM /dH . The peak effect in the magnetization 
is not reported at H'^^^ , but is observed well below H^^ .^^^ 

Magnetostriction measurements have also been per- 
formed as a function of H.^^ The magnetostriction 
changes the slope at Hp indicating a second order tran- 
sition. Magnetostriction has also been measured in an 
oblique field. The FFLO regime shrinks in an oblique 
field and when H is tilted by 20° from the ah plane, the 
FFLO regime occupies only a very narrow regime below 
the upper critical field. 

5.3.3 Thermal conductivity 

Thermal conductivity is a unique transport quan- 
tity that does not vanish in the superconducting state. 
It is particularly suitable for probing the quasiparti- 
cle structure, responding only to dclocalized low energy 
quasiparticle excitations, since localized quasiparticles 
do not carry heat. The thermal conductivity is a direc- 
tional probe of quasiparticles. In recent years, the ther- 
mal conductivity has been used effectively to probe the 
anisotropy of the order parameter in unconventional su- 
perconductors."* This is due to the fact that normal dclo- 
calized quasiparticles are easily excited along the nodal 
directions. In the FFLO state, anisotropy of the dclo- 
calized quasiparticle excitation spectrum is particularly 
noticeable, because the electron velocity vanishes along 
certain directions. The inherent anisotropy of the quasi- 
particle excitation is expected to produce anisotropy of 
the thermal conductivity. 

Thermal conductivity measurements in a parallel field 
have been carried out in two different configurations, H \\ 
j and H _L j, near i?j2-^^ Here j is the heat current. The 
thermal conductivity as a function of 77 for Jf || j clearly 
displays a kink at H^^ and exhibits an enhancement above 
In contrast, it is difficult to resolve these features for 

The observed enhancement of the thermal conductiv- 
ity for iJ II J is contrary to initial expectations, since in 
this geometry the FFLO nodal planes are created per- 
pendicular to the heat current, and therefore enhance- 
ment of the thermal conductivity would not be expected. 
A possible scenario is the formation of a peculiar quasi- 
particle structure within the vortex core in CeCoIns, as 
will be discussed in §6.2. Then, as suggested in Ref. [22], 
the vortex core size might increase at the nodal planes. 
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Fig. 10. Transverse sound velocity for two different configu- 
rations, the Lorentz mode, || fc || [100], u ||[010], and non- 
Lorentz mode, || fJ" || [100], k ||[010]. In the non-Lorentz mode, 
the sound wave couples only to the crystal lattice, while the 
Lorentz mode corresponds to a flux line bending mode. The dif- 
ference between the two modes can be regarded as the response 
from the flux line lattice. 



leading to enhancement of the thermal conductivity. Fur- 
ther theoretical investigations into the thermal transport 
of the inhomogcneous FFLO state are necessary. 

5.3.4 Ultrasound velocity 

In order to study possible FFLO states, it is particu- 
larly important to determine the structure of the fiux 
line lattice (FLL), which is intimately related to the 
electronic structure. At the FFLO transition, the elas- 
tic modulus of the vortex lattice is expected to change, 
which can be detected by ultrasound measurements. In 
ultrasound experiments, sound waves are coupled to the 
FLL when it is pinned by crystal lattice defects. The 
ensuing modified sound dispersion allows detailed infor- 
mation about the FLL-crystal coupling to be extracted. 
More precisely, the change in the sound velocity with 
respect to the sound velocity in the normal state is^*^ 

A«t = , . „^!/o,,,o, . , (19) 



fioPmVt 1 + {^oBy/^\k\^/fiojc' 

where pm is the mass density of the material, k is the 
wave vector of the sound, and jc is the critical current 
density. 

Figure 10 displays the transverse sound velocity for 
H II a6 as a function of temperature for two different 
configurations with respect to the polarization vector m, 
the sound propagation vector fc, and H.^° We distin- 
guish between the Lorentz force mode (Fig. 10, u _L H), 
where the sound wave couples to the vortices through 
the induced Lorentz force Fl ~ X~^{u x H) x B, and 
the non-Lorentz mode {u \\ H), where the sound wave 
couples only to the crystal lattice. The transverse sound 
velocity in the Lorentz mode, in which the flux motion 
is parallel to the a6-plane, is strongly enhanced com- 
pared to that in the non-Lorentz mode. This indicates 
that the transverse ultrasound strongly couples to the 
FLL in the Lorentz mode. The difference between the 
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sound velocities in the Lorentz and non-Lorentz modes, 
Avt = (tt -L H) — v^{u II H), can be regarded as be- 
ing the contribution of the flux line lattice to the sound 
velocity, which provides direct information of C44 of the 
FLL. 

The enhancement of in the non-Lorentz mode in- 
dicates a stiffening of the crystal lattice at the super- 
conducting transition below T^^H). The enhancement of 
the sound velocity in the Lorentz mode is much greater 
than that in the non-Lorentz mode. This immediately 
indicates that as the temperature is lowered the FLL 
becomes pinned. In the Lorentz mode, the ultrasound 
velocity changes its slope with a cusp at T*^ . It should be 
emphasized that such a cusp structure is not observed 
in the non-Lorentz mode. This provides direct evidence 
that the transition at Tj* is characterized by a change in 
flux line pinning. The sound velocity below Tj* is smaller 
than that extrapolated from above T^^^ indicating that 
pinning is weaker below the FFLO transition. 

Since the volume pinning force = j^B is determined 
by the balance of the energy gain obtained by placing a 
flux line on a crystal defect, and the energy loss due to 
the ensuing deformation of the FLL, a decrease of pin- 
ning has two possible origins. Namely, the decrease of the 
elementary pinning interaction between the defects and 
the flux lines, and the increase of one or several of the 
FLL elastic moduli. In the present case, the relevant elas- 
tic modulus is the dispersive flux line tilt modulus C44. 
While a full theory is as yet unavailable, it is expected 
that the coupling of the deformations of the FLL and the 
nodal planes will lead to a slight decrease of C44 in the 
FFLO phase. ^"'^^ This would, however, lead to stronger 
pinning, as the FLL would be able to better adapt to 
the pinning impurities. The relatively smaller sound ve- 
locity in the FFLO state suggests that the first effect, 
a decrease of the elementary pinning interaction due to 
depression of the average order parameter magnitude, as 
indicated by the penetration depth and magnetization 
measurements, is the main effect. 

Thus ultrasound velocity measurements reveal an im- 
portant decrease of the superconducting volume fraction 
at the FFLO transition. The nature and magnitude of 
the ultrasound velocity anomaly is quantitatively consis- 
tent with the modulated order parameter predicted for 
the FFLO phase.20'118 

5.3.5 Nuclear magnetic resonance 

One of the most powerful probes for extracting micro- 
scopic information on quasiparticle and spin excitations 
is NMR. Here we discuss the NMR results in accordance 
with Ref. [24], in which a dramatic change of the NMR 
spectra at {T*^,H*t^) has been reported. 

Figure 11 displays the NMR spectra at the In(l) site 
(^^^In (/ = 9/2)) with axial symmetry in the Celna layer, 
which is located in the center of the square lattice of Ce 
atoms. (The tetragonal structure of CeCoIns consists of 
alternating layers of Celna and CoIn2 and hence has two 
inequivalent In sites per unit cell.) The ^^^In NMR spec- 
trum is plotted as a function of Knight shift "^^^K, which 
is a local spin susceptibility at the In(l)-site, obtained 
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Fig. 11. Main panel: ^-"^^In-NMR spectra outside, slightly above 
H^i (blue line), slightly above T* (black), and well inside (red) 
the FFLO phase. The resonance feature at higher frequency is 
marked by hatching, (i?, T)-points at which each NMR spec- 
trum was measured are shown by crosses in the inset. Inset: 
Experimental H-T phase diagram for CeCoIns below 1 K for 
a II ah. The regions shown in yellow and blue depict the FFLO 
and BCS regimes, respectively. Horizontal arrows indicate the 
magnetic fields at which the NMR spectra were measured (see 
also Fig. 12(a)). 



from the central ^^^In line (±1/2 ^ Tl/2 transition). 
For comparison, the spectrum in the normal state slightly 
above (blue fine), in the BCS state slightly above tJ' 
(black line), and well inside the FFLO state (red line), 
are shown in Fig. 11. 

A remarkable feature is the appearance of a new reso- 
nance peak with small but finite intensity at higher fre- 
quency in the spectrum well inside the FFLO phase, as 
clearly seen at ^^^K ~ 2.03% (hatched region). It is no- 
table that the Knight shift of this new higher resonance 
peak coincides well with that of the resonance line in the 
normal state above B^.^. (blue), while the position of the 
lower resonance line is located close to that of the BCS 
state above T|*(7J) (black). Therefore, it is natural to 
deduce that the higher resonance line originates from a 
normal quasiparticle regime, which forms below T^J{H). 
The appearance of the higher new resonance line, which 
is a manifestation of a normal quasiparticle structure, is 
what is exactly expected for the FFLO state. It is re- 
ported that the intensity of the higher resonance line is 
about 3-5 percent of the total intensity and is nearly 
T-independent below 180 mK, indicating that the total 
volume of the planar node occupies a few percent of the 
total volume.^"* 

The well separation of the two resonance lines defi- 
nitely indicates that the quasiparticles excited around 
the FFLO planar nodes are spatially well-separated from 
those excited around vortex cores. This is again consis- 
tent with the topological structure of the vortex in the 
FFLO state, as discussed in §2.5, i.e., the quasiparticle 
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Fig. 12. ^^^In-NMR spectra as a function of Knight shift from 
the normal state (black lines) to the superconducting state (red 
lines) at (a) H = 11.3 T for H \\ ab (see the inset of Fig. 11) and 
(b) H = 4.80 T for H 1. ab. The intensity is normalized by the 
largest peak intensity. 



regions within the vortex hne do not spatiaUy overlap 
with those in the FFLO planar nodes. 

Figure 12(a) displays the temperature evolution of the 
spectra at H = 11.3 T when crossing the {T^ , H*^)-lme 
(see the horizontal arrow in the inset of Fig. 11). The 
higher resonance line grows rapidly with decreasing T 
just below T^*{H). A double peak structure appears at 
T = 240 mK, followed by a shoulder structure at T = 
260 ~ 300 mK, indicating that the intensity of the higher 
resonance line dominates. The two lines merge into a 
single line above T*{H) ~ 300 mK. 

It has been suggested that the temperature evolu- 
tion of the NMR spectra is compatible with the quasi- 
particle nodal planes formed in the FFLO state. 
The field induced layered structure expected in the 
FFLO phase resembles the superconducting states of 
stacks of superconductor-normal-superconductor (S-N- 
S) Josephson tunnel junctions. In NMR experi- 
ments, rf shielding supercurrents flow across the FFLO 
planar nodes. Because of the second order transition 
at T^{H), the modulation length of the order param- 
eter, i.e., the thickness of the superconducting layers, 
A (= 2T:/q), diverges as A cx (T* - T)"" with a > 
upon approaching T^*{H). Therefore A exceeds the in- 
plane penetration length in the vicinity of T|*. In this 
situation, an rf field penetrates into the normal sheets 
much deeper than into the SC sheets, which results in a 
strong enhancement of the NMR intensity from the nor- 
mal sheets. At low temperature where A becomes com- 
parable to ^ (<g; A) , an rf field penetrates the normal 
sheets to the same degree as the superconducting sheets. 
The temperature evolution of the higher resonance line 
can be explained semi-quantitatively by this model. 

To date, three groups have performed NMR ex- 
periments to investigate the FFLO state in parallel 
fields.^''' Although all groups have reported changes 
of the spectra in the vicinity of i?|2i the features of 



the spectra are different for each group and the inter- 
pretations of the spectra are controversial. Ref. [28] re- 
ports a change of the line width of the main spectrum 
at (Tipi?!*) without the appearance of the higher reso- 
nance peak. In Ref. [120], NMR spectra at the In(l) and 
In(2) sites were measured in a wide frequency range. Al- 
though no discernible change of the spectra at the In(l) 
site below is observed, the appearance of a broad line 
at the In(2) site has been reported. This broadening is 
attributed to possible antiferromagnetic (AFM) ordering 
within the vortex cores, which will be discussed in §5.6.2. 
At present, the origin of the difference in the NMR spec- 
tra among the three groups is an open question and calls 
for further investigation. 

5.4 FFLO state in a perpendicular field 

Despite the observation of an anomaly in the heat ca- 
pacity in a perpendicular field just below H^, the study 
of the superconducting state has focused on the H \\ ab 
case. Recent NMR and magnetization measurements pro- 
vide several pieces of evidence for the FFLO state in the 
perpendicular field case as well. NMR results in a perpen- 
dicular field have been reported by only one group. Fig- 
ure 12(b) depicts the temperature evolution of the NMR 
spectrum for H ± ab just below (4.95 T at T = 0). 
The spectrum and its temperature evolution bear strik- 
ing resemblance to those for H \\ ab shown in Fig. 12(b). 
In the superconducting state, the spectrum exhibits a 
shoulder structure at lower frequencies (T — 0.46 K). 
This shoulder structure is followed by a double peak 
structure as the temperature decreases, clearly showing 
below T = 0.43 K. The double peak structure persists 
down to the lowest temperature (70 mK) and the inten- 
sity of the higher resonance line decreases with decreas- 
ing temperature. As discussed in §5.3.5, the appearance 
of the higher resonance peak is a manifestation of a new 
quasiparticle region in the FFLO state. Therefore, the 
NMR data for H J- ab also leads to the conclusion that 
the FFLO state is realized for H _L ab. 

According to Ref. [25], the double peak anomaly in 
the spectrum is not observed below H = 4.7 T, indicat- 
ing that the critical field FF^ separating the FFLO and 
the BCS pairing state lies in a very narrow range be- 
tween 4.7 T and 4.75 T. Very recent magnetization mea- 
surements reports that dM/dH exhibits a peak at H'^. 
The results indicate that the phase boundary between 
FFLO and BCS states is nearly horizontal, i.e., i?^ is F- 
independent, which is consistent with the NMR results 
(see Fig. 9(a)). 

Very recently, NMR spectra at the Co, In(l), and In(2) 
sites have been measured in a wide frequency range for 
H _L ab.^^^ In contrast to the spectra for H \\ ab re- 
ported in Ref. [120], in which a possible AFM order- 
ing is suggested within the vortex cores as discussed in 
§5.3.5, no anomalous broadening of the spectra has been 
observed in this configuration. Therefore, there is no mi- 
crosopic evidence of the AFM ordering in the vortex core 
for H _L ab. 
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5.5 H-T phase diagram 

We here discuss the FFLO region in the experimentally 
determined H-T phase diagram displayed in Figs. 9(a) 
and 9(b). The topology of the phase transition lines for 
both configurations H \\ ab and H ^- ah is identical 
in that the FFLO line branches from the first order to 
normal transition line. Despite the similarity of the phase 
diagrams of both configurations, there are distinguishing 
differences between them. For H ± afe, the FFLO phase 
occupies a tiny high-iJ/low-T corner in the H-T phase 
diagram, in contrast to the case for H \\ ab. In addition, 
H*^ has almost no temperature dependence, in striking 
contrast to ilf || ah. 

The appearance of the FFLO state in H A- ab is not 
surprising, because as discussed in §5.1, the Maki pa- 
rameter a'^ =5.0 greatly exceeds the minimum value of 
1.8 for the FFLO state. We note that is even larger 
than a". The extremely large can be attributed to 
the strongly enhanced Pauli paramagnetic susceptibility 
in the normal state. In fact, the Pauli susceptibility for 
H ^- ah \s nearly twice as large as that for H \\ ah. 

The sizable decrease of the FFLO region for H J. ah 
indicates that the FFLO state is more stable for || ah. 
This may be related to the quasi-2D Fermi surface of 
the system, as discussed in §2.3. Recently, it has been 
proposed that the presence of a small sheet of 3D Fermi 
surface can stabilize the FFLO state in a perpendicu- 
lar field. In addition, the enhanced antiferromagnetic 
fluctuation in a perpendicular field compared with a par- 
allel field may be important for the stability of the FFLO 
state. 

Some of FFLO theories assume a second order tran- 
sition at the upper critical field, which is incompati- 
ble with the first order transition observed in CcCoIns. 
The first order transition from the normal state to the 
FFLO state has been theoretically predicted by several 

authors.38.40.47,56.75,76,100,123 _ 

Recently, the phase diagram in a quasi-2D system with 
rf-wave symmetry has been examined in the context of 
CeCoIn5.^5,i24 rpj^g second order FFLO-BCS transition 
that branches from the first order transition line at the 
upper critical field is reproduced in Ref. [75] , as shown in 
Fig. 6. However, an FFLO state with Landau level index 
n = 1 at very low temperatures has not been reported 
experimentally. Recently it has been pointed out that the 
FFLO state with n = 1 is more unstable against the im- 
purities than that with n ~ 0.^^^ For the second-order 
transition, the i/Jj-curve has been semi-quantitatively 
reproduced. In addition, it has been pointed out that 
the fluctuation-driven first order transition by impuri- 
ties is consistent with observations for CeCoIns.^^' 
The effect of the Fermi liquid correction has also been 
discussed in connection with CeCoIn5.^°° 

5.6 Another aspect of the FFLO state in CeCoIn^ 

5. 6. 1 Influence of antiferromagnetic fl,uctuation 

The above results strongly suggest that an FFLO state 
exists in CeCoIns for both parallel and perpendicular 
fields. However, several unexpected features have been 
observed. A most interesting and important feature is 



the competition or coexistence with AFM ordering. 

The occurrence of superconductivity in CeCoIns in the 
vicinity of the AFM quantum critical point (QCP) gives 
rise to pronounced non-Fermi-liquid behavior in various 
physical quantities in the normal state. 127-129 More- 
over, the specific heat and resistivity measurements have 
shown a crossover from Fermi liquid to non-Fermi liquid 
behavior in a magnetic field, which indicates the pres- 
ence of a field- tuned AFM QCP at or in the vicinity of 
H^.^^''^'^''''^ This implies that superconductivity persists, 
preventing the development of AFM order. Moreover, it 
has been reported that AFM ordering appears when a 
magnetic field is applied in the superconducting state of 
CeRhlns, which exhibits superconductivity under pres- 
sure. 

Then an important question is whether the origin 
of the the second order transitions at {T^,H*^) and 
{T^jH"^) are related to some magnetic transitions. Ex- 
perimentally, the coexistence of an FFLO state and an 
AFM ordering in CeCoIns has been a controversial issue, 
as discussed below. 

For H II ah, as discussed in §5.3.5, recent measure- 
ments of NMR spectra report the appearance of a re- 
markable broad line at the In(2) site, which has been ten- 
tatively interpreted as an AFM ordering within the vor- 
tex cores, -'^2° as observed in high-Tc cuprates.^^^' ^'^^ The 
coexistance of an AFM ordering and an FFLO state is 
appeared to be possible, since an FFLO q-vector is much 
smaller than an AFM wave number. On the other hand, 
there are several pieces of evidence that is against the 
AFM ordering in the high field superconducting phase 
of CeCoIns. First, recent heat capacity measurements 
of CeCoIns under pressure clearly demonstrate that the 
FFLO regime expands in the H-T plane when the AFM 
fluctuations are suppressed by pressure. This provides di- 
rect evidence that the AFM fluctuations are unfavorable 
for the stability of the FFLO state. This is understood by 
considering that the FFLO state favors forming the para- 
magnetic regime, while the AFM fluctuations that prefer 
the singlet state prevent the formation of the paramag- 
netic regime. Second, under AFM order, the alternating 
hyperfine fields would produce two inequivalent ^^^In(l) 
sites, which would give rise to two NMR resonance lines 
with equal intensities, but such resonance lines are not 
observed, as shown in Figs. 11, 12(a), and 12(b). Third, 
an AFM transition generally leads to a large anomaly 
of the stiffness constant of the crystal lattice, which can 
be detected by ultrasound. However, as shown in Fig. 
10, no anomaly of the ultrasound is observed at T|* in 
the no-Lorentz mode. Fourth, no signature of the tran- 
sition line has been observed above -ffj2i unlike the case 
of CeRhlns. These considerations imply that the AFM 
ordered moments induced by magnetic field in the vortex 
state for H \\ ah are very tiny, if present. 

For H _L ah, no anomalous broadening of the NMR 
spectra has been observed as discussed in §5.4,1^^ in con- 
trast to if _L ah. In addition, the nearly T-indcpendent 
horizontal FFLO-BCS boundary in the H-T plane for 
H J- ah shown in Fig. 9(a) is very different from the AFM 
order line reported in Ref. [131]. Therefore, at least for 
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H ± o,h, there is no evidence of the AFM ordering in the 
vortex state. 

Until now, there is no theory of what might be ex- 
pected for the FFLO state if magnetism wore induced 
by a magnetic field. Further experimental and theoreti- 
cal investigation is strongly required to clarify this issue. 

5.6.2 Unusual vortex core 

It has been suggested that the NMR spectra also pro- 
vide an important information on the nature of the vor- 
tex core, which in turn would be closely related to the 
FFLO state. The main NMR spectra just below H^^2 
(black and red lines in Fig. 11) deviates significantly from 
that in the normal state above H^}2- Since this regime is 
located at ~ 1/3 of the orbital limited upper critical 
field, a substantial portion (~ 1/3) of the NMR intensity 
should originate from the vortex core. Therefore a dis- 
tinct deviation of the NMR spectra just below ijjg from 
the normal state indicates that the Knight shift within 
the vortex cores is significantly reduced from the normal 
value. This feature is in sharp contrast to that in con- 
ventional superconductors, where the Knight shift within 
cores coincides with that in the normal state. 

As discussed in §5.3.5 and §5.6.1, a possible AFM 
vortex core has been suggested, as observed in high- 
Tc cuprates.^^^'^^^ This peculiar feature of the vortex 
core in CeCoIns is also inferred from the energy scale of 
the quasiparticlo spectrum in the vortex core set by the 
confinement energy hcjo ~ /^a^^^ 88, 134-136 almost 
all superconductors, this energy level is of the order of 
fjK — mK, which is negligibly small. For CeCoIn5, how- 
ever, ef ^ 15 K and A ^ 5K, so that hojQ ^ 1.5 K 
and the vortex core spectrum becomes important at low 
temperatures, i.e., the system is in a quantum regime. 
The small ep is due to the rcnormalized Fermi energy by 
the heavy electron mass. Moreover, in CeCoIns, £ ~ 3 /um 
and ^ ~ 4 nm, which yields £/^ ~ 750 at low fields, which 
is much larger than ep/A ^ 3. This is in sharp contrast 
to other superconductors in which £/^ <C £f/A. This im- 
plies that the vortex in CeCoIn5 enters a new quantum 
regime, the superclcan regime, that is difficult to access 
in most superconductors.^* 

Although the possibility of achieving the superclean 
regime and AFM ordering within cores in CeCoIns 
should be scrutinized, both effects are expected to af- 
fect the quasiparticle structure and phase diagram in the 
FFLO state. Further detailed investigations are neces- 
sary to clarify this issue. 

6. Summary 

Several heavy fermion and layered organic supercon- 
ductors are considered good candidates for exhibiting 
the FFLO state. Particularly, the peak effect observed 
in the magnetization near Hc2 has been attributed to 
the FFLO state in several heavy fermion superconduc- 
tors. However, subsequent studies have raised questions 
about the interpretation of the data. Although a pos- 
sible FFLO state has been suggested in several organic 
superconductors in parallel fields, no thermodynamic or 
microscopic evidence has been reported so far. Of these 



superconductors, CeCoIus appears to be a unique system 
for which corroborating experimental evidence supports 
the FFLO state. In fact, CeCoIn5 appears to meet all the 
strict requirements on the existence of the FFLO state. 
The FFLO state in CeCoIn5 is present for both parallel 
and perpendicular field directions. The topology of the 
phase transitions in H-T phase diagram for both config- 
urations appear to be fully determined. However, some 
of the microscopic measurements are still controversial 
and some of the results seem to be different from those 
originally predicted for the FFLO state in some respects. 

For further clarification of the FFLO state, exper- 
iments using more direct techniques, such as neutron 
diffraction that can detect the spatial variation of the 
polarized spin excitations, and scanning tunneling mi- 
croscopy that can directly detect the spatial variation of 
the quasiparticle density of states, are required. More- 
over phase sensitive experiments, such as measuring the 
tunnel effect between the FFLO state and the conven- 
tional s-wave BCS state, would provide an unambigu- 
ous method to detect the FFLO state. For instance, it 
has been suggested that while the Josephson current be- 
tween FFLO and conventional BCS superconductors is 
suppressed, it can be recovered by applying a magnetic 
field. It has also been suggested that tunneling spectra 
between FFLO and conventional BCS superconductors 
have characteristic line shapes with several maxima and 
minima reflecting minigap structures due to the periodic 
pair potentials. ^•^^ Recently, it has been pointed out that 
the intrinsic pinning of vortices in the FFLO state can 
be used as a direct probe of the spatial modulation of 
the FFLO state. 

Although not discussed in this paper, interesting states 
which resemble the FFLO state have been proposed for 
various systems. For example, an FFLO-like state is in- 
duced near the interface between a ferromagnet and a 
superconductor. In this case, the pair amplitude in the 
ferromagnetic region spatially oscillates and as a result, a 
O-TT transition is expected. (For a review, see Ref. [140].) 
In superconductors with no spatial inversion symmetry, 
two novel superconducting states are predicted, depend- 
ing on the strength of the spin-orbit coupling Aso- When 
Aso ^ A, a helical and a stripe vortex states have been 
proposed, in which the order parameters possesse a mod- 
ulation e"'"^^^^ and cos(q' • r),^^^ respectively, which are 
similar to the FFLO state. When Ago ^ A, a possible 
triplet FFLO state induced by spin-orbit coupling in 
the absence of an external magnetic field has been pro- 
posed. 

Acknowledgments 

We thank C. Capan, R. Ikeda, U. Klein, K. Kuma- 
gai, K. Machida, K. Maki, K. Nagai, D. Rainer, and 
C. J. van der Beek for many clarifying discussions on 
the subject of this paper. We also thank H. Adachi, 
T. Hanaguri, M. Ichioka, K. Izawa, K. Kakuyanagi, 
Y. Kasahara, R. Movshovich, M. Nohara, Y. Onuki, 
T. Sakakibara, T. Shibauchi, S. Shishido, M. Tachiki, 
M. Tanatar, Y. Tanaka, S. Uji, and T. Watanabe for 
helpful discussions. This work was partly supported by a 
Grant-in- Aid for Scientific Reserch from the Ministry of 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 15 



Education, Culture, Sports, Science and Technology. 



M. Sigrist and K. Ucda: Rev. Mod. Phys. 63 (1991) 239. 
V. P. Mineev and K. V. Samokhin; in Introduction to Un- 
conventional Superconductivity (Gordon and Breach Science 
Publishers, 1999). 

P. Thalmeicr and G. Zwicknagl: in Handbook on the Physics 
and Chemistry of Rare Earths, Vol.34, (North Holland 
2005). 

For a recent review, see Y. Matsuda, K. Izawa, and 
I. Vekhter: J. Phys.; Condens. Matter 18 (2006) R705. 
P. Fuldc and R. A. Ferrell: Phys. Rev. 135 (1964) A550. 
A. I. Larkin and Yu. N. Ovchinnikov: Zh. Eksp. Tcor. Fiz. 47 
(1964) 1136 [translation: Sov. Phys. JETP 20 (1965) 762]. 
For a recent review, see R. Casalbuoni and G. Narduli; Rev. 
Mod. Phys. 76 (2004) 263. 

L. W. Gruenberg and L. Gunther; Phys. Rev. Lett. 16 (1966) 
996. 

L. G. Aslamazov: Sov. Phys. JETP 28 (1969) 773. 

S. Takada: Prog. Theor. Phys. 43 (1970) 27. 

K. Gloos, R. Modler, H. Schimanski, C. D. Brcdl, C. Geibel, 

F. Steglich, A. I. Buzdin, N. Sato, and T. Komatsubara: 
Phys. Rev. Lett. 70 (1993) 501. 

G. Yin and K. Maki: Phys. Rev. B 48 (1993) 650. 

H. Burkhardt and D. Rainer: Ann. Physik 3 (1994) 181. 
H. Shimahara; Phys. Rev. B 50 (1994) 12760. 

N. Dupuis: Phys. Rev. B 51 (1995) 9074. 

M. Taxihiki, S. Takahashi, P. Gegenwart, M. Weiden, 

M. Lang, C. Geibel, F. Steglich, R. Modler, C. Paulsen, and 

Y. Onuki: Z. Phys. B 100 (1996) 369. 

Rad03 

A. I. Buzdin and H. Kachkachi: Phys. Lett. A 225 (1997) 

341. 

H. A. Radovan, N. A. Fortune, T. P. Murphy, S. T. Hannahs, 
E. C. Palm, S. W. Tozer, and D. Hall: Nature 425 (2003) 51. 
A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and 
J. L. Sarrao: Phys. Rev. Lett. 91 (2003) 187004. 
T. Watanabe, Y. Kasahara, K. Izawa, T. Sakakibara, Y. Mat- 
suda, C. .1. van der Beck, T. Hanaguri, H. Shishido, R. Settai, 
and Y. Onuki: Phys. Rev. B 70 (2004) 020506(R). 
T. Watanabe, K. Izawa, Y. Kasahara, Y. Haga, Y. Onuki, 
P. Thalmeicr, K. Maki, and Y. Matsuda: Phys. Rev. B 70 
(2004) 184502. 

C. Capan, A. Bianchi, R. Movshovich, A. D. Christianson, 

A. Malinowski, M. F. Hundley, A. Lacerda, P. G. Pagliuso, 
and J. L. Sarrao: Phys. Rev. B 70 (2004) 134513. 

V. F. Correa, T. P. Murphy, C. Martin, K. M. Pur- 
cell, E. C. Palm, G. M. Schmiedeshoff, J. C. Cooley, and 
S. W. Tozer: cond-mat/0609487. 

K. Kakuyanagi, M. Saitoh, K. Kumagai, S. Takashima, 
M. Nohara, H. Takagi, and Y. Matsuda: Phys. Rev. Lett: 
94 (2005) 047602. 

K. Kumagai, M. Saitoh, T. Oyaizu, Y. Furukawa, 
S. Takashima, M. Nohara, H. Takagi, and Y. Matsuda: to 
be published in Phys. Rev. Lett.; cond- mat/0605394. 

C. F. Miclea, M. Nicklas, D. Parker, K. Maki, J. L. Sarrao, 
J. D. Thompson, G. Sparn, and F. Steglich: Phys. Rev. Lett. 
96 (2006) 117001. 

X. Gratcns, L. M. Ferrcira, Y. Kopelcvich, N. F. Olivcira Jr., 
P. G. Pagliuso, R. Movshovich, R. R. Urbano, J. L. Sarrao, 
and J. D. Thompson: cond-mat/0608722. 
V. F. Mitrovic, M. Horvatic, C. Berthier, G. Knebel, 
G. Lapertot, and J. Flouquet: cond-mat/0605305. 

D. Saint-James, G. Sarma, and E. J. Thomas: Type II Su- 
perconductivity (Pergamon, New York, 1969). 

B. S. Chandrasekhar, Apply. Phys. Lett. 1 (1962) 7; 
A.M. Clogston, Phys. Rev. Lett. 9 (1962) 266. 

J.B. Kettcrson and S.N. Song: Superconductivity (Cambridge 
University Press, 1999). 

S. Takada and T. Izuyama: Prog. Theor. Phys. 41 (1969) 635. 
L. N. Bulaevskii: Zh. Eksp. Teor. Fiz. 65 (1973) 1278; [trans- 



34 

35 

36; 
37 

38 

39; 

40 
41 

42 

44' 
45 

46; 
47; 

48; 

49; 

50 

51 

52 

53 
54 

55 

56 

57; 
58' 



59 
60 
61 
62 
63 
64 
65 
66; 
67 

68 

69; 
70 
71 

72 

73 
74 



75 

76; 
77; 

78) 



lation: Sov. Phys. JETP 38 (1974) 634]. 

K. Aoi, W. Dieterich, and P. Fulde: Z. Physik 267 (1974) 

223. 

Y. Suzumura and K. Ishino: Prog. Theor. Phys. 70 (1983) 
654. 

K. Machida and H. Nakanishi: Phys. Rev. B 30 (1984) 122. 

H. Shimahara: J. Phys. Soc. Jpn. 67 (1998) 736. 

J. A. Bowers and K. Rajagopal: Phys. Rev. D 66 (2002) 

065002. 

C. Mora and R. Combescot: Europhys. Lett. 66 (2004) 833; 
R. Combescot and G. Tonini: Phys. Rev. B 72 (2005) 094513. 

C. Mora and R. Combescot: Phys. Rev. B 71 (2005) 214504. 
Q. Wang, H.-Y. Chen, C.-R. Hu, and C. S. Ting: Phys. Rev. 
Lett. 96 (2006) 117006. 

R. Combescot and G. Tonini: Phys. Rev. B 72 (2005) 094513. 
K. Maki and H. Won: Czech. J. Phys. 46 (1996) 1035. 
K. Yang and S. L. Sondhi: Phys. Rev. B 57 (1998) 8566. 
A. B. Vorontsov, J. A. Sauls, and M. J. Graf: Phys. Rev. B 
72 (2005) 184501. 

R. Ikeda and H. Adachi: Phys. Rev. B 69 (2004) 212506. 

D. F. Agtcrbcrg and K. Yang: J. Phys.: Condens. Matter 13 
(2001) 9259. 

C. Caroli, P. G. de Gennes, and J. Matricon, Phys. Lett. 9 
(1964) 307. 

T. Mizushima, K. Machida, and M. Ichioka: Phys. Rev. Lett. 

95 (2005) 117003. 

T. Mizushima, K. Machida, and M. Ichioka: Phys. Rev. Lett. 
94 (2005) 060404. 

H. Takayama, Y. R. Lin-Liu, and K. Maki: Phys. Rev. B 21 
(1980) 2388. 

H. Shimahara: J. Phys. Soc. Jpn. 67 (1998) 1872; Physica B 
259-261 (1999) 492. 

Y. Ghashi: J. Phys. Soc. Jpn. 71 (2002) 2625. 

The QLRO is sustained in the states expressed by cqs. (6), 

(7), and (8). See Rcf. [52]. 

K. V. Samokhin and M. S. Mar'cnko: Phys. Rev. B 73 (2006) 
144502. 

D. Dalidovichi and K. Yang: Phys. Rev. Lett. 93 (2004) 
247002. 

R. Modler et al: Phys. Rev. Lett. 76 (1996) 1292. 

A. Yamashita, K. Ishii, T. Yokoo, J. Akimitsu, M. Hedo, 

Y. Inada. Y. Onuki, E. Yamamoto, Y. Haga, and R. Kadono: 

Phys. Rev. Lett. 79 (1997) 3771. 

T. Kiss et al: Phys. Rev. Lett. 94 (2005) 057001. 

T. Sakakibara et al.: a preprint. 

C. Geibel et al: Z. Phys. B 84 (1991) 1. 

N. K. Sato et al: Nature 410 (2001) 340. 

M. Jourdan, M. Huth, and H. Adrian: Nature 398 (1999) 47. 

N. Bernhoeft: Europhys. J. B 13 (2000) 685. 

J. R. Clem: Phys. Rev. B 43 (1991) 7837. 

M. R. Norman: Phys. Rev. Lett. 71 (1993) 3391. 

V. Vinokur, B. Khaykovich, E. Zeldov, M. Konczykowski, 

R. A. Doyle, and P. H. Kes: Physica C 295 (1998) 209. 

M. B. GaifuUin, Y. Matsuda, N. Chikumoto, J. Shimoyama, 

and K. Kishio: Phys. Rev. Lett. 84 (2000) 2945. 

N. Avraham et al: Nature 411 (2001) 451. 

S. S. Banerjee et al: Phys. Rev. B 58 (1998) 995. 

A. M. Troyanovski: J. Aarts, P. H. Kes: Nature 399 (1999) 

665. 

H. R. Ott, H. Rudigier, T. M. Rice, K. Ueda, Z. Fisk, and 
J. L. Smith: Phys. Rev. Lett. 52 (1984) 1915. 
H. Tou et al.: unpublished. 

H. Shimahara, S. Matsuo, and K. Nagai: Phys. Rev. B 53 
(1996) 12284: Y. Suginishi and H. Shimahara: Phys. Rev. B 

74 (2006) 024518. 

H. Adachi and R. Ikeda: Phys. Rev. B 68 (2003) 184510. 

M. Houzet and V. P. Mineev: cond-mat/0606799. 

H. Shimahara and D. Rainer: J. Phys. Soc. Jpn. 66 (1997) 

3591. 

A. I. Buzdin and J. P. Brison: Europhys. Lett. 35 (1996) 
707; Phys. Lett. A 218 (1996) 359. A gap function of the 
form A(r) oc exp [— imip] cxp [— p^/2]/9'" with integer m is 
assumed, m is not the Landau level index n. 



16 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



79) U. Klein, D. Rainer, and H. Shimahaxa: J. Low Temp. Phys. 
118 (2000) 91. 

80) M. Houzet and A. Buzdin: Europhys. Lett. 50 (2000) 375. 

81) U. Klein: Phys. Rev. B 69 (2004) 134518. 

82) K. Yang and A. H. MacDonald: Phys. Rev. B 70 (2004) 
094512. 

83) C. Petrovic, R. Movshovich, M. Jaime, P. G. Pagliuso, 
M. F. Hundley, J. L. Sarrao, Z. Fisk, and J. D. Thompson: 
Europhys. Lett. 53 (2001) 354. 

84) R. Settai, H. Shishido, S. Ikeda, Y. Murakawa, 
M. Nakashima, D. Aoki, Y. Haga, H. Harima, and Y. 
Onuki: J. Phys.: Condens. Matter 13 (2001) L627. 

85) Y. Kawasaki, S. Kawasaki, M. Yashima, T. Mito, 

G. -Q. Zheng, Y. Kitaoka, H. Shishido, R. Settai, Y. Haga, 
and Y. Onuki: J. Phys. Soc. Jpn. 72 (2003) 2308. 

86) R. Movshovich, M. Jaime, J. D. Thompson, C. Petrovic, 
Z. Fisk, P. G. Pagliuso, and J. L. Sarrao: Phys. Rev. Lett. 
86 (2001) 5152. 

87) K. Izawa, H. Yarnaguclii, Y. Matsuda, H. Shishido, R. Settai, 
and Y. Onuki: Phys. Rev. Lett. 87 (2001) 057002. 

88) Y. Kasahara, Y. Nakajima, K. Izawa, Y. Matsuda, 
K. Behnia, H. Shishido, R. Settai, and Y. Onuki: Phys. Rev. 
B 72 (2005) 214515. 

89) N. J. Curro, B. Simovic, P. C. Hammel, P. G. Pagliuso, 
J. L. Sarrao, J. D. Thompson, and G. B. Martins: Phys. 
Rev. B 64 (2001) 180514. 

90) Y. Kohori, Y. Yamato, Y. Iwamoto, T. Kohara, E. D. Bauer, 
M. B. Maple, and J. L. Sarrao: Phys. Rev. B 64 (2001) 
134526. 

91) S. Ozcan, D. M. Broun, B. Morgan, R. K. W. Haselwimmer, 
J. R. Waldram, J. L. Sarrao, S. Kamal, C. P. Bidinosti, and 
P. J. Turner: Europhys. Lett. 62 (2003) 412; R. J. Ormeno, 
A. Sibley, C. E. Gough, S. Sebastian, and L R. Fisher: 
Phys. Rev. Lett. 88 (2002) 047005; E. E. M. Chia, D. J. Van 
Harlingen, M. B. Salamon, B. D. Yanoff, I. Bonalde, and 
J. L. Sarrao: Phys. Rev. B 67 (2003) 014527. 

92) P. M. C. Romke, M. A. Tanatar, C. S. Turel, J. Berdeklis, 
C. Petrovic, and J. Y. T. Wei: Phys. Rev. Lett. 94 (2005) 
107005; W. K. Park, L. H. Greene, J. L. Sarrao, and 
J. D. Thompson: cond-mat/0606535. 

93) H. Aoki, T. Sakakibara, H. Shishido, R. Settai, Y. Onuki, 
P. Miranovic, and K. Machida: J. Phys.: Condens. Matter 
16 (2004) L13. 

94) A. Vorontsov and L Vekhter: Phys. Rev. Lett. 96 (2006) 
237001. 

95) T. Tayama, A. Harita, T. Sakakibara, Y. Haga, H. Shishido, 
R. Settai, and Y. Onuki: Phys. Rev. B 65 (2002) 180504. 

96) C. Martin, C. C. Agosta, S. W. Tozer, H. A. Radovan, 

E. C. Palm, T. P. Murphy, and J. L. Sarrao: Phys. Rev. 
B 71 (2005) 020503(R). 

97) N. Oeschler, P. Gegenwart, M. Lang, R. Movshovich, 
J. L. Sarrao, J. D. Thompson, and F. Steglich: Phys. Rev. 
Lett. 91 (2003) 076402. 

98) A. Bianchi, R. Movshovich, N. Oeschler, P. Gegenwart, 

F. Steglich, J. D. Thompson, P. G. Pagliuso, and 
J. L. Sarrao: Phys. Rev. Lett. 89 (2002) 137002. 

99) L. DeBeer-Sclimitt, C. D. Dewhurst, B. W. Hoogenboom, 
C. Petrovic, and M. R. Eskildsen: Phys. Rev. Lett. 97 (2006) 
127001. 

100) A. B. Vorontsov and M. J. Graf: cond-mat/0606572. 

101) The Pauli limiting field H^2 reduced by a factor \/xo/x> 
while 11^2^'-' is reduced by Xo/x- 

102) H. Shimahara: J. Phys. Soc. Jpn. 66 (1997) 541; J. Phys. Soc. 
Jpn. 68 (1999) 3069. 

103) T. Ishiguro: J. Phys. IV France 10 (2000) Pr3-139. 

104) M. A. Tanatar, T. Ishiguro, H. Tanaka, and H. Kobayashi: 
Phys. Rev. B 66 (2002) 134503. 

105) S. Uji, H. Shinagawa, T. Terashima, T. Yakabe, Y. Terai, 
M. Tokumoto, A. Kobayashi, H. TanaJja, and H. Kobayashi: 
Nature 410 (2001) 908. 

106) L. Balicas, J. S. Brooks, K. Storr, S. Uji, M. Tokumoto, 

H. TanaJja, H. Kobayashi, A. Kobayashi, V. Barzykin, and 
L. P. Gor'kov: Phys. Rev. Lett. 87 (2001) 067002. 



107) H. Shimahara: J. Phys. Soc. Jpn. 71 (2002) 1644; Physica B 
329-333 (2003) 1442. 

108) M. Houzet, A. Buzdin, L. Bulaevskii, and M. Maley: Phys. 
Rev. Lett. 88 (2002) 227001. 

109) S. Manalo and U. Klein: J. Phys.: Condens. Matter 12 (2000) 
L471. 

110) S. Uji, T. Terashima, M. Nishimura, Y. Takahide, 
T. Konoike, K. Enomoto, H. Cui, H. Kobayashi, 
A. Kobayashi, H. Tanaka, M. Tokumoto, E. S. Choi, 
T. Tokumoto, D. Graf, and J. S. Brooks: Phys. Rev. Lett. 97 
(2006) 157001. 

111) M.-S. Nam, J. A. Symington, J. Singleton, S. J. Blundell, 
A. Ardavan, J. A. A. J. Perenboom, M. Kurmoo, and P. Day: 
J. Phys.: Condens. Matter 11 (1999) L477-L484. 

112) J. Singleton, J. A. Symington, M.-S. Nam, A. Ardavan, 
M. Kurmoo, and P. Day: J. Phys.: Condens. Matter 12 (2000) 
L641; J. A. Symington, J. Singleton, M.-S. Nam, A. Arda- 
van, M. Kurmoo, and P. Day: Physica B 294-295 (2001) 
418. 

113) K. Izawa, H. Yamaguclii, T. Sasaki, and Y. Matsuda: Phys. 
Rev. B 88 (2002) 027002. 

114) Physics of Organic Superconductors and Conductors (edited 
by A. G. Lebed, Springer, Berlin, to be published in 2007). 

115) M. Decroux and O. Fisher: in Superconductivity in Ternary 
Compounds. Part II (Springer, Heidelberg, 1982). 

116) M. Houzet and A. Buzdin: Phys. Rev. B 63 (2001) 184521. 

117) D. Dominguez, L. Bulaevskii, B. Ivlev, M. Maley, and 
A. R. Bishop: Phys. Rev. B 53 (1996) 6682. 

118) R. Ikeda: cond-mat/0610796. 

119) Y. Matsuda, M. B. Gaifullin, K. Kumagai, K. Kadowaki, and 
T. Mochiku: Phys. Rev. Lett. 75 (1995) 4512. 

120) B.-L. Young, R. R. Urbano, N. J. Curro, J. D. Thompson, 
J. L. Saxrax3, A. B. Vorontsov, and M. J. Graf: cond- 
mat/0608040. 

121) T. Oyaizu, K. Kumagai, H. Shishido, and Y. Matsuda: un- 
published. 

122) R. Ikeda: eond-mat/0610863. 

123) S. Matsuo, S. Higashitani, Y. Nagato, and K. Nagai: J. Phys. 
Soc. Jpn. 67 (1998) 280. 

124) H. Won, K. Maki, S. Haas, N. Oeschler, F. Weickert, and 
P. Gegenwart: Phys. Rev. B 69 (2004) 180504(R). 

125) R. Ikeda and H. Adachi: Phys. Rev. Lett. 95 (2005) 269703. 

126) D. Dalidovichi and K. Yang: Phys. Rev. Lett. 95 (2005) 
269704. 

127) A. Bianchi, R. Movshovich, I. Vekhter, P. G. Pagliuso, and 
J. L. Sarrao: Phys. Rev. Lett: 91 (2003) 257001. 

128) Y. Nakajima et al: J. Phys. Soc. Jpn. 73 (2004) 5; ibid 75 
(2006) 023705. 

129) V. A. Sidorov, M. Nicklas, P. G. Pagliuso, J. L. Sarrao, 
Y. Bang, A. V. Balatsky and J. D. Thompson: Phys. Rev. 
Lett. 89 (2002) 157004. 

130) J. Paglione, M. A. Tanatar, D. G. Hawthorn, E. Boaknin, 
R. W. Hill, F. Ronning, M. Sutherland, L. Taillefer, C. Petro- 
vic, and P. C. Canfield: Phys. Rev. Lett. 91 (2003) 246405. 

131) T. Park, F. Ronning, H. Q. Yuan, M. B. Salamon, 
R. Movshovich, J. R. Sarrao, and J. D. Thompson: Nature 
440 (2006) 65. 

132) K. Kakuyanagi, K. Kumagai, Y. Matsuda, and M. Hasegawa: 
Phys. Rev. Lett. 90 (2003) 197003. 

133) B. Lake et al: Science 291 (2001) 1759. 

134) Y. Matsuda, N. P. Ong, Y. F. Yan, J. M. Harris, and J. B. Pe- 
terson: Phys. Rev. B 49 (1994) 4380. 

135) J. M. Harris et al: Phys. Rev. Lett. 73 (1994) 1711. 

136) K. Kumagai, K. Nozaki, and Y. Matsuda: Phys. Rev. B 63 
(2001) 144502. 

137) K. Yang and D. F. Agterberg: Phys. Rev. Lett. 84 (2000) 
4970. 

138) Y. Tanaka, Y. Asano, M. Ichioka, and S. Kashiwaya: cond- 
mat/0607665. 

139) L. Bulaevskii, A. Buzdin, and M. Maley: Phys. Rev. Lett. 90 
(2003) 067003. 

140) A. I. Buzdin: Rev. Mod. Phys. 77 (2005) 935. 

141) R. P. Kaur, D. F. Agterberg, and M. Sigrist: Phys. Rev. Lett. 



J. Phys. Soc. Jpn. Full Paper Author Name 17 

94 (2005) 137002. 143) H. Tanaka, H. Kaneyasu, and Y. Hasegawa: cond- 

142) D. F. Agterberg, and R. P. Kaur: cond-mat/0612216. mat/0609751. 



